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ABSTRACT 

A  real-space  molecular-orbital  description  of  electronic  wavefunctions 
which  are  postulated  to  be  the  precursors  of  the  superconducting  state  in 
metals,  alloys,  compounds,  and  noncrystalline  materials  Is  presented,  based 
on  self-consistent-field  X-alpha  scattered-wave  (SCF-Xa-SW)  molecular-orbital 
calculations  for  clusters  representing  the  local  molecular  environments  in 
these  materials.  It  is  shown  through  a  variety  of  examples  that  there  is  a 
persistent  correlation  between  the  occurrence  of  superconductivity  in  a 
material  and  the  existence  of  molecular  orbitals  at  the  Fermi  energy  with 
coherent  (i.e.,  in-phase)  spatially  extended  (e.g.,  pir,  da,  dz2,  or  so) 
bonding  atomic-orbital  components.  It  is  also  argued  that,  while  these 
individual  molecular-orbital  components  are  usually  "one-dimensional"  or 
"two-dimensional"  in  nature,  the  composite  precursor  superconducting  state 
of  a  material  arises  from  a  coherent  three-dimensional  "network"  or  "array" 
of  these  components.  It  is  further  demonstrated  through  examples  that  for 
each  such  orbital  component,  covalent  bond  overlap  in  a  line  or  plane  of 
atoms  is  enhanced  by  the  interaction  of  atoms  not  in  that  line  or  plane,  as 
shown  by  the  pir-po*  bonding-antibonding  interaction  in  aluminum  and  by 
metal-ligand  interactions  in  A15  and  layered  transition-metal  compounds. 

This  description  of  the  precursor  superconducting  state  is  consistent  with 
the  original  conjectures  of  London  and  Slater  that  the  superconducting-state 
wavefunction  is  "molecular"  in  nature,  "rigid"  in  character,  and  of  wide 
spatial  extent,  from  which  observed  physical  properties  (e.g.,  diamagnetism 
and  nondissipative  electrical  currents)  of  the  superconducting  state  logically 
follow.  This  molecular-orbital  description  is  further  argued  to  be  consistent 
with  Cooper's  concept  of  electron  pairing  in  the  superconducting  state  through 
a  net  attractive  electron-electron  interaction  but  differs  from  the  usual 
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viewpoint  In  attributing  this  attractive  Interaction  directly  to  real -space 
chemical -bond-1  Ike  formation  rather  than  to  electron-phonon  coupling  per  se. 
The  two  points  of  view,  however,  are  not  necessarily  Incompatible.  The 
above  molecular-orbital  criteria  for  superconductivity  are  diametrical  to 
those  for  the  occurrence  of  local  magnetic  moments  and  ferro-  and  anti  ferro¬ 
magnetism,  namely  the  existence  of  spin-polarized  spatially  localized 
anti  bonding  (e.g.,  do*  and  dir*)  molecular  orbitals  near  the  Fermi  energy. 
These  criteria  therefore  provide  a  simple  conceptual  basis  for  understanding 
the  generally  mutually  exclusive  incidence  of  superconductivity  and  magnetism 
among  the  elements  of  the  periodic  table,  although  they  can  also  be  used  to 
explain  the  occasional  coexistence  of  superconductivity  and  ferromagnetism 
or  anti  ferromagnetism  in  some  materials.  The  molecular-orbital  criteria  for 
superconductivity  complement  existing  formal  theories  in  that  they  appear 
to  explain  simply  and  directly  certain  observed  chemical  trends  of  super¬ 
conductivity,  such  as  those  frequently  emphasized  by  Matthias,  and  are 
naturally  applicable  to  superconducting  materials  lacking  long-range  crystal¬ 
line  order,  such  as  solid-solution  alloys,  amorphous  alloys,  and  small 
particles.  Finally,  it  appears  that  the  molecular-orbital  approach  can  be 
used  to  explain  from  simple  considerations  why  some  materials  (e.g.,  Cu,  Ag, 
and  Au)  are  neither  superconducting  nor  magnetic,  to  rule  out  certain  types 
of  substances  (e.g.,  certain  classes  of  quasi- one-dimensional  organic  solids; 
metallic  hydrogen  at  attainable  high  pressures)  as  possible  superconductors, 
and  to  suggest  ways  of  systematically  improving  existing  classes  or 
synthesizing  new  classes  of  superconducting  materials. 
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I .  INTRODUCTION 

The  Bardeen-Cooper-Schrieffer  (BCS)  theory  of  superconductivity^ 

ascribes  the  onset  of  the  superconducting  state  of  a  crystal  at  the  transition 

temperature,  Tc>  to  electrons  attractively  paired  via  thermal  vibrations  of 

the  lattice  (phonons).  The  BCS  formula  for  T£  is 

Tc*  eD  exp  [-1/VNUf)],  (1) 

where  eD  is  the  Debye  temperature,  V  is  the  net  attractive  potential  between 

2 

electrons  in  Cooper  pairs  induced  by  the  electron-phonon  interaction,  and 
N( e p)  is  the  electronic  density  of  states  at  the  Fermi  energy.  The  BCS 
theory  together  with  its  strong-coupling  extension  has  been  eminently  suc¬ 
cessful  in  accounting  for  the  physical  properties,  e.g.,  nondissipative 
current,  diamagnetism,  and  thermodynamics  of  the  superconducting  state,  and 
has  correlated  many  experimental  data  in  terms  of  a  few  basic  parameters. 

3 

Nevertheless,  it  has  often  been  emphasized  by  Matthias  that  the  BCS  theory 
and  formal  extensions  thereof  do  not  satisfactorily  explain  the  observed 
dependence  of  superconductivity  on  crystal  structure  and  chemistry 
(especially  for  transition  metals,  alloys,  and  compounds),  do  not  address 
the  mutually  exclusive  incidence  of  superconductivity  and  magnetism  for  the 
elements  in  the  periodic  table,  and  are  not  very  useful  for  predicting  new 
superconducting  materials.  The  existence  of  superconducting  materials 
possessing  only  short-range  structural  order,  such  as  superconducting 

A  ° 

amorphous  alloys  and  superconducting  metal  particles  down  to  *  50  A  in 
5  6 

size,  also,  strictly  speaking,  lies  outside  the  domain  of  BCS  theory. 

The  BCS  theory,  in  its  conventional  form,  is  based  on  the  traditional 
concepts  of  long-range  crystalline  order  and  momentum  (£)  space,  leading  to 
coherence  lengths  of  the  superconducting  state  that  greatly  exceed  the 
short-range  order  and  electron  mean  free  path  characteristic  of  amorphous 
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and  small -particle  superconductors.  This  suggests  the  need  for  a  local 

chemical -bonding  or  "real-space"  molecular  description  of  the  superconducting 

7  8 

state.  Indeed,  London  *  In  his  phenomenological  approach  to  superconductivity 
discusses  the  possibility  of  developing  a  molecular  description  of  the  super- 
conducting  state  (see  Chapter  E  of  Ref.  8),  and  Slater  in  his  classic  papers 
on  the  subject  provides  a  qualitative  description  of  the  spatial  character 
of  the  superconducting-state  wavefunctlon.  With  speculations  that  mechanisms 
other  than  electron -phonon  coupling  can  attractively  pair  electrons  in  the 
superconducting  state, a  molecular  criterion  that  accounts  for  the 
known  chemical  trends  in  the  occurrence  of  superconductivity  would  be  a 
useful  tool  in  the  continuing  effort  to  identify  new  classes  of  superconducting 
materials. 

3 

The  systematic  trends  of  Tc  with  valence-electron-per-atom  ratio  and 
the  difficulty  of  calculating  BCS  parameters  for  transition  metals,  their 
alloys,  and  compounds  have  prompted  the  development  of  rudimentary  chemical 
models,^2-14  based  on  atomic  d  orbitals  and  aimed  primarily  at  the  estima¬ 
tion  of  Tc.  However,  these  models  are  not  easily  generalized  to  non¬ 
transition-metal  superconductors  and  do  not  yield  a  real-space  molecular 

description  of  the  superconducting  state.  The  most  general  chemical  approach 

1 5 

to  superconductivity  is  that  of  Krebs,  who  has  suggested  that  superconduc¬ 
tivity  is  possible  only  if  the  normal  chemical  bonding  system  in  the  crystal 
or  parts  of  the  crystal  permits  the  construction  of  a  molecular  wavefunction 
which  for  at  least  one  space  direction  is  not  intersected  by  plane  or  conical 
nodal  surfaces,  and  if  the  corresponding  electron  band  is  not  fully  occupied. 
Note  that  Krebs  attempts  only  to  account  for  the  incidence  of  superconductiv¬ 
ity  through  qualitative  chemical  bonding  arguments  and  not  the  evaluation  of 
a  specific  electron-pairing  mechanism  and  T  . 

In  this  paper,  we  present  a  real-space  molecular-orbital  description 


-6- 


and  criteria  for  the  occurrence  of  the  superconducting  state  in  metals, 

alloys,  compounds,  and  noncrystalline  materials,  including  amorphous  solids 

and  small  particles.  This  theoretical  approach  to  superconductivity  makes 

direct  use  of  the  results  of  first-principles  self-consistent-field  X-alpha 

scattered-wave  (SCF-Xa-SW)  molecular-orbital  calculations  for  clusters 

representing  the  local  molecular  environments  in  these  materials.  The 

SCF-Xa-SW  method1®  is  based  on  the  combined  use  of  the  Xa  density-functional 

approximation1 7  to  electron-electron  exchange  and  correlation  and  the  multiple- 

1 8 

scattered-wave  technique  of  solving  the  Schrodinger  self-consistent- field 

molecular-orbital  equations.  In  conjunction  with  the  "transition-state" 

procedure1®’17  for  calculating  electronic  transitions  between  orbitals, 

this  method  leads  to  a  description  of  both  the  ground  and  excited  electronic 

states  of  a  system.  The  calculation  of  different  molecular  orbitals  for 

different  spins  allows  one  to  consider  the  effects  of  magnetic  spin  polariza- 

19  20 

tion.  Along  with  applications  to  a  wide  variety  of  polyatomic  molecules,  * 

21-23  24 

including  biological  macromolecules  and  polymers,  this  approach  has 

been  used  successfully  to  elucidate  the  local  bulk  and  surface  electronic 

25-27 

structures  and  properties  of  solids,  e.g.,  crystalline  semiconductors 

28-31  32 

and  localized  defects  and  impurities  therein,  amorphous  semiconductors, 

33-38 

crystalline  and  amorphous  metals  and  alloys  including  magnetic 

35  37  24  39  40 

properties,  ’  molecular  crystals,  and  chemisorption  on  metal  surfaces  ’ 

41  42 

including  the  calculation  of  photoelectron  emission  intensities,  * 

43  44 

electron-scattering  cross  sections,  and  "shake-up"  spectra  for  chemi¬ 
sorbed  species.  These  theoretical  studies  collectively  suggest  that  short- 
range  order  and  local  chemical  bonding  largely  determine  the  measurable 
electronic  structures  and  properties  of  many  types  of  materials,  whether 
ordered  or  disordered. 

In  view  of  the  complementary  natures  of  superconductivity  and  magnetism, 
it  should  be  noted  that  spin-unrestricted  SCF-Xa-SW  cluster  molecular-orbital 
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studies  Indicate  that  the  local  magnetic  moments  of  Fe  and  Mn  Impurities 

in  crystalline  copper  above  the  Kondo  temperature  are  associated  with  the 

existence  of  spin-polarized  anti  bonding  (a*  and  **)  Fe(d)-Cu(s,d)  molecular 

orbitals  near  the  Fermi  energy  (see  Figs.  10,  12,  and  13  of  Ref.  35).  Similar 

17  26 

cluster  models  for  ferromagnetic  bcc  crystalline  o-lron  *  suggest  that  the 

onset  of  collective  magnetism  below  the  Curie  temperature  and  the  persistence 

of  "spin  clusters”  above  the  Curie  temperature  are  associated  with  the 

existence  of  spin-polarized  antibonding  (a*  and  it*)  Fe(d)-Fe(d)  molecular 

orbitals  near  the  Fermi  energy  of  the  type  shown  in  Fig.  1.  Recent  spin- 

polarized  SCF-Xa-SW  cluster  models  for  amorphous  iron-boron  and  i ron-phosphorus 
38 

alloys  are  also  in  good  quantitative  agreement  with  measured  magnetic  properties 
and  clearly  show  the  key  role  of  anti  bonding  molecular  orbitals  at  the  Fermi 
energy  in  determining  how  alloying  alters  the  magnetic  properties  of  iron. 


For  the  sake  of  completeness,  it  may  also  be  noted  that  the  SCF-Xa-SW  method 

22 

has  been  applied  to  iron-containing  biological  molecules  such  as  ferredoxin 
23 

and  hemoglobin  and  the  results  used  to  explain  observed  magnetic  properties 
(including  the  calculation  of  magnetic  hyperfine  parameters  )  as  functions 
of  temperature  and  oxygenation,  respectively.  Finally,  the  spin-unrestricted 
SCF-Xa-SW  method,  in  conjunction  with  the  transition-state  procedure J 6,1 ^ 
has  recently  been  employed  with  good  results  to  calculate,  from  first  principles, 
ferromagnetic  and  anti ferromagneti c  exchange  coupling  constants  (Heisenberg 

i  45 

exchange  integrals)  for  transition-metal  cluster  coordination  complexes. 

The  most  general  finding  of  the  theoretical  studies  of  magnetic  substances 
summarized  above  is  that  both  local  ami  collective  forms  of  magnetism  are 
associated  with  the  occurrence  of  anti bonding  spin-polarized  molecular 
orbitals  (e.g.,  the  da*  and  dir*  components  in  Fig.  1)  near  the  Fermi  energy. 
Because  of  the  complementary  natures  of  magnetism  and  superconductivity, 

:j 


t; 
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it  is  reasonable  to  expect  that  SCF-Xa-SW  cluster  molecular-orbital  studies 

of  superconducting  materials  can  also  provide  insight  into  the  role  of 

chemical  bonding  in  determining  the  superconducting  state.  Indeed,  it  will 

be  shown  in  this  paper  that  the  molecular-orbital  criteria  for  the  occurrence 

of  superconductivity  and  magnetism  are  diametrical,  thereby  offering 

a  simple  explanation  for  the  mutually  exclusive  occurrence  of  these  phenomena 

in  a  material  at  the  same  temperature  and  pressure.  The  orbital  criteria  for 

2 

superconductivity  will  further  be  shown  to  be  consistent  with  Cooper's 
concept  of  electron  pairing  through  a  net  attractive  electron-electron 
interaction  (V  in  Eq.  1)  but  differs  from  BCS  theory^  in  attributing  this 
attractive  interaction  directly  to  real-space  chemical-bond  formation  at 
the  Fermi  energy,  rather  than  to  electron -phonon  coupling  per  se.  Finally 
the  molecular-orbital  approach  presented  herein  will  be  used  to  explain 
why  some  metals  (e.g. ,  Cu,  Ag,  and  Au)  are  neither  superconducting  nor  magnetic 
and  to  provide  a  reliable  theoretical  tool  for  predicting  new  superconducting 
materials . 

II.  CLUSTER  MOLECULAR  ORBITALS  AND  BAND  STRUCTURE 

The  relationships  of  SCF-Xa-SW  cluster  molecular-orbital  energy  levels 
and  densities  of  states  to  the  corresponding  crystal  band  structures  and 
densities  of  states  calculated  from  Bloch's  theorem  and  measured  spectro¬ 
scopically  have  been  discussed  in  previous  publications  and  reports  for  the 
examples  of  copper,  ’  3  nickel,  aluminum,  ’  y  and  iron.  Nevertheless, 
in  order  to  clarify  the  use  of  cluster  molecular  orbitals  as  a  basis  for 
superconductivity,  we  point  out  some  specific  similarities  and  differences 
between  cluster  electronic  structure  and  crystal  band  structure. 

For  simplicity,  consider  first  the  example  of  an  isolated  benzene 
molecule,  the  pit  orbitals  of  which  are  shown  schematically  on  the  left  side 
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of  Fig.  2.  The  tt^  orbital,  which  is  occupied  by  two  spin-paired  electrons 
and  has  the  energy  level  ,  is  the  most  strongly  bonding  of  the  pir  molecular 
orbitals,  whereas  the  degenerate  ^  and  tt^  orbitals,  which  together  are 
occupied  by  four  spin-paired  electrons  and  have  the  energy  level  3*  are 
relatively  weakly  bonding.  The  unoccupied  orbital  is  strongly  antibonding, 
while  the  unoccupied  and  orbitals  are  weakly  antibonding.  Because 
of  the  one-dimensional  periodicity  of  the  carbon  atoms  around  the  benzei-r 
ring,  the  molecular  orbitals  can  be  expressed  in  symmetrized  Bloch  for;;., *8 
resulting  in  the  discrete  one-dimensional  "k-space"  band  profile  shown 
schematically  on  the  right  side  of  Fig.  2.  The  "Fermi  energy"  ep  separates 
the  occupied  states  from  the  unoccupied  ones.  In  this  particular  example, 
there  is  a  simple  one-to-one  correspondence  between  the  "real-space"  molecular 
orbitals  and  "k-space"  band  structure,  as  indicated  by  the  connecting 
lines  in  Fig.  2.  It  has  often  been  noted  that  the  diamagnetic  "ring  currents" 
induced  in  benzene  and  other  aromatic  molecules  through  the  response  of  the 
pir  electrons  to  an  external  magnetic  field  are  nondissipative  currents  similar 
in  some  (but  not  all)  respects  to  the  persistent  currents  of  a  superconducting 
ring.8’47*48 


On  the  other  hand,  the  molecular  orbitals  of  a  finite  cluster  repre¬ 
senting  the  local  molecular  environment  in  a  crystal  transform  according  to 
the  irreducible  representations  of  the  cluster  point  symmetry  group  and 
therefore  are  not  generally  equivalent  on  a  one-to-one  basis  to  Bloch  band 
eigenstates  of  the  crystal,  except  at  the  points  in  Ic-space  which  have  the 
full  cluster  point  group  symmetry.  For  example,  the  F(l<=o)  Bloch  state  has 
the  full  0^  point-group  symmetry  in  a  cubic  crystal.  Nevertheless,  it  is 
possible  to  construct  combinations  of  degenerate  Bloch  states  of  different 
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ic  (the  so-called  "star"  of  the  wavevector)  that  transform  according  to  the 

49 

irreducible  representations  of  the  cluster  point  group.  From  this,  one 
can  conclude  that  there  should  be  a  reasonably  close  correspondence  between 
cluster  and  crystal  densities  of  states.  Examples  of  such  correspondence  will 
be  found  in  Refs.  33,  34,  37,  and  39. 

As  an  illustrative  application  of  the  above  concepts  to  a  superconducting 

metal,  consider  the  example  of  aluminum.  In  Fig.  3,  the  SCF-Xct-SW  molecular- 

orbital  energy  levels  of  a  19-atom  cluster  representing  the  local  environment 

34 

of  fee  crystalline  aluminum  up  to  second-nearest  neighbors  are  compared  with 
the  energy  bands  of  the  crystal  shown  along  principal  symmetry  directions 
of  the  Brillouin  zone.  The  Al-jg  cluster  density  of  states,  based  on 
a  Gaussian  broadening  of  the  discrete  molecular-orbital  energy  levels  of  the 
type  described  in  Refs.  33  and  34,  is  shown  in  Fig.  4.  A  detailed  comparison 
of  cluster  and  crystal  densities  of  states  for  aluminum,  including  a  cluster 
molecular-orbital  analysis  of  the  density  of  states  derived  f v om  photoemission 
spectra,  is  presented  in  Ref.  34. 

Like  the  results  described  in  Ref.  35  for  a  19-atom  cluster  representing 
crystalline  copper,  there  is  a  remarkably  close  correspondence  of  the  total 
energy  span  of  the  Al-jg  cluster  orbitals  to  the  band  width  of  the  crystal 
and  a  close  relationship  between  individual  cluster  orbitals  to  particular  energy- 
band  Bloch  eigenstates  or  combinations  thereof.  For  example,  the  lowest 
occupied  cluster  energy  level,  la^,  shown  in  Fig.  3  is  the  discrete  molecular- 
orbital  analogue  of  the  I^(i?=o)  s-band  Bloch  state.  Among  the  highest  occupied 
and  lowest  unoccupied  cluster  orbitals,  the  fully  occupied  tlg  level  corresponds 
mainly  to  p-band  states  just  below  the  Fermi  energy  e^,  the  partially 
occupied  ^t-^Up)  level  corresponds  to  combinations  of  degenerate  and 

A-j  p-band  states  at  the  Fermi  energy,  and  the  empty  t2g  level  corresponds 
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state.  58  Similar  studies  of  Cr  for  Interatomic  distances  corresponding  to 


primarily  to  I ^  p-band  states  just  above  the  Fermi  energy.  The  tj  ,  4t^ u( e p ) > 

and  2tgg  molecular  orbitals  are  responsible  for  the  peak  in  the  cluster 

density  of  states  at  the  Fermi  energy  evident  in  Fig.  4  and  are  the  discrete 

analogues  of  the  parallel  and  z-|  p  bands  largely  responsible  for  the  peak 

in  the  crystalline  band-structure  and  photoemission  densities  of  states  around 
34 

the  Fermi  energy. 

As  will  be  shown  in  the  following  section  of  the  paper,  the  t1u(ep) 
cluster  molecular  orbital  is  key  to  understanding  the  real-space  nature  of 
the  superconducting  state  of  aluminum.  The  wavefunction  contour  maps  of 
this  orbital,  plotted  in  the  (200)  and  (110)  crystallographic  planes,  are 
shown  in  Figs.  5  and  6,  respectively.  The  solid  and  dashed  contours  represent 
positive  and  negative  regions,  respectively,  of  the  molecular-orbital  wave- 
function.  Regions  of  net  overlap  between  contours  of  the  same  sign  are 
bonding,  whereas  regions  of  positive  and  negative  contours  separated  by  nodes 
are  anti  bonding.  The  most  striking  feature  of  Figs.  5  and  6  is  the  "one¬ 
dimensional"  ir-bonding  character  of  the  A1  p  orbitals  along  the  "horizontal" 
direction  of  each  figure,  much  like  the  one-dimensional  p  orbitals  of  the 
benzene  molecule  shown  in  Fig.  2.  However,  in  contrast  to  benzene,  where 
the  pir i  molecular  orbital  lies  well  below  the  Fermi  energy,  the  pir-bonding 
orbital  of  aluminum  is  located  exactly  at  the  Fermi  energy  at  T  =  0°K. 

Unlike  benzene,  moreover,  there  is  also  a  "one-dimensional"  a*-antibonding 
component  between  the  A1  p  orbitals  along  the  "vertical"  direction  in  Figs.  5 
and  6,  which,  through  electron-electron  repulsion,  isolates  the  pir  orbitals 
along  parallel  lines  of  atoms  and  synergistically  promotes  covalent  pn  bond 
overlap  within  each  line  of  atoms.  Although  Figs.  5  and  6  imply  anisotropy  in 
each  t-| u(£ p)  molecular-orbital  component,  the  triple  degeneracy  of  this 
orbital  leads  to  equivalent  pir-bonding  and  pa *-anti bonding  components  along 
the  two  complementary  orthogonal  crystal  directions  and  therefore  to  an 
isotropic  three-dimensional  orbital  charge  distribution.  As  cluster  size 
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is  systematically  increased  to  include  more  shells  of  atoms  in  the  crystal 
(e.g.,  see  the  43-atom  results  in  Ref.  34),  each  pir-bonding  and  po*- 
anti bonding  orbital  component  at  the  Fermi  energy  expands  spatially  and,  in 
the  limit  of  a  large  number  of  atoms,  forms  a  coherent  macroscopic  molecular 
orbital  of  the  type  shown  in  Fig.  7.  A  perspective  drawing  of  pir  bonding 
along  a  line  of  atoms  is  shown  in  Fig.  8.  This  “real-space"  character  of 
aluminum  molecular  orbitals  at  the  Fermi  energy  cannot  be  extracted  easily 
from  the  conventional  "Ic-space"  Bloch  description  of  the  crystal  band 
structure,  except  in  hindsight  as  an  appropriate  combination  of  E-| ,  e3,  and 
Bloch  eigenstates,  and  thus  has  eluded  solid-state  physics. 

III.  CLUSTER  MOLECULAR  ORBITALS  AND  THE  SUPERCONDUCTING  STATE 

The  traditional  view  of  the  BCS  and  earlier  theories  of  superconducti vity^  ’ 

is  that  the  superconducting  state  of  a  crystal  is  a  "special,"  highly  ordered 

electronic  state  of  large  but  finite  coherence  length,  constructed  from  linear 

combinations  of  degenerate  Bloch  states  of  opposite  momenta  or  wavevectors 

(t,  -£)  at  the  Fermi  energy  and  occupied  by  spin-paired  electrons  between 

which  there  is  a  net  attraction.  In  this  section  of  the  paper,  we  will 

attempt  to  establish  that  the  above  view  is  consistent  with  the  real-space 

representation  of  the  superconducting  state  as  arising  from  discrete  molecular 

orbitals  at  the  Fermi  energy  with  coherent,  spatially  extended  bonding 

components,  for  which  certain  cluster  molecular  orbitals  of  the  type  described 

in  the  preceding  section  (the  ptr-pa*  orbitals  of  aluminum)  are  precursors. 

One  attractive  feature  of  this  concept  is  that  it  leads  immediately  to  a 

simple  explanation  of  the  aforementioned  similarity  between  the  nondissipative 

diamagnetic  ring  currents  associated  with  the  pir  electrons  of  benzene  (see 

8  47  48 

Fig.  2)  and  other  aromatic  molecules  *  *  and  the  persistent  diamagnetic 

currents  in  a  superconductor.  Moreover,  since  the  cluster  molecular-orbital 
model  is  not  dependent  on  the  assumption  of  long-range  crystalline  order, 
it  is  applicable  to  amorphous^  and  small -parti cle^’^  superconductors. 
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where  "£  is  not  a  good  quantum  number,  strictly  speaking. 

In  the  above  picture,  the  attractive  interaction  responsible  for  electron 

pairing  in  the  superconducting  state  occurs  through  the  exclusive  occupancy 

at  the  Fermi  energy  (for  T  ^  0°K)  of  spatially  extended  covalent  bonding 

molecular-orbital  components  of  the  p-n-  type  for  p-electron  materials  like 

aluminum  or,  as  will  be  shown  below,  bonding  orbitals  of  the  d<5  and  dz2  type 

for  transition  metals  and  their  compounds.  This  is  not  to  imply  that  the 

phonons  do  not  play  a  role  in  superconductivity,  since  the  phonon  spectrum  and 

the  concomitant  electron-phonon  interaction  are  byproducts  of  the  attractive 

potential  energy  of  the  atoms  associated  with  the  chemical  bond.  The 

electron  pairs  are  spatially  delocalized  over  the  "coherence  length"  of  the 

appropriate  bonding  molecular  orbital ,  e.g. ,  the  pir  components  in  Figs.  5,  6, 

and  7,  in  much  the  same  fash’on  as  the  electron  pairs  are  delocalized  over 

the  benzene  pn-j -electron  ring  in  Fig.  2.  Moreover,  the  electron-electron 

repulsion  associated  with  the  pc*-anti bonding  components  of  Figs.  5  through  7 

synergistically  promotes  the  covalent  pairing  of  electrons  associated  with 

2 

the  pir-bonding  components.  Note  that  in  the  original  paper  by  Cooper  on 

the  concept  of  electron  pairing  in  superconductivity  and  in  the  much  earlier 
g 

papers  by  Slater  on  the  nature  of  the  superconducting  state,  the  attractive 

interaction  responsible  for  the  pairing  of  electrons  of  opposite  Bloch 

wavevectors  {t,  -t)  at  the  Fermi  energy  can,  in  retrospect,  be  identified 

directly  with  chemical  bonding,  without  altering  the  basic  ideas  and 

respective  formalisms  in  these  papers, 
g 

Slater  has  also  shown  that  if  the  superconducting  state  is  characterized 
by  a  discrete  energy  level  with  a  reasonably  uniform  wavefunction  which  is 
spatially  delocalized  over  ^  137  atomic  diameters  and  then  decays  exponentially 
in  a  molecular  fashion,  this  state  acts  like  a  huge  "atom"  or  "molecule"  and 
exhibits,  for  moderate  applied  magnetic  fields,  a  large  diamagnetism  of  the 
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7  8 

type  observed  in  superconductors  and  predicted  by  London's  *  phenomenological 

theory  of  superconductivity.  The  cluster  ptr-pc*  molecular-orbital  wavefunction 

of  aluminum  described  above  is  very  much  like  the  wavefunction  described  by 
g 

Slater,  especially  as  the  cluster  size  is  expanded  to  encompass  many  atoms 
(see  Fig.  7).  The  "wide  extension  in  space"  of  the  superconducting-state 
wavefunction  is  also  discussed  by  London  (see  p.  150  of  Ref.  8).  The  electron 
pairs  occupying  this  molecular  orbital  should  not  be  viewed  as  propagating 
free  electrons  in  the  sense  of  ordinary  conducting  electrons,  but  rather  as 
rigidly  confined  to  the  coherence  length  of  the  wavefunction  (a  hundred  or 
more  atoms,  but  not  the  entire  crystal),  like  electrons  confined  as  "standing 
deBroglie  waves"  in  a  large  box  of  constant  potential  energy.  The  wide  spatial 
extent  of  the  superconducting-state  molecular  orbitals  suggests  from  the 
uncertainty  principle  that  the  electron  pairs  occupying  this  orbital  exist 

O 

in  the  same  center-of-mass  momentum  state.  As  discussed  by  London  and 

g 

Slater,  such  a  wavefunction  is  unchanged  by  the  application  of  an  external 
magnetic  field  of  moderate  size.  In  other  words,  the  electrons  in  this 
state,  which  represent  a  small  fraction  of  the  total  number  of  electrons  in 
the  system,  are  prevented  from  "coiling  up"  like  free  electrons  when  the 
system  is  brought  into  a  magnetic  field,  but  stay  essentially  as  they  are 

Q 

without  the  magnetic  field,  as  if  "frozen."  The  synergistic  pit-bonding  and 
po*- anti  bonding  components  of  the  molecular-orbital  wavefunction  shown  in 
Fig.  7  are  responsible  for  the  "rigid"  nature  of  the  superconducting  state 
of  aluminum. 

It  is  appropriate  here  to  point  out  an  interesting  analogy  between 
the  present  real -space  concept  of  superconductivity  and  the  coherent  trans¬ 
mission  of  light  by  optical  fibers.  The  spatially  extended  pit-bonding 
molecular-orbital  components  of  aluminum  shown  in  Fig.  7  (and,  as  will  be 
shown  below,  the  dfi-bonding  orbitals  of  transition  metals)  can  be  viewed 


simplistically  as  providing,  through  covalent  bond  overlap,  quasi-one- 
dimensional  "waveguides"  for  the  standing  deBroglie  waves  associated  with 
the  coherent  superconducting-state  electrons,  in  analogy  to  the  way  optical 
fibers  act  as  electromagnetic  waveguides  for  the  coherent  transmission  of 
light.  Indeed,  the  orthogonal  pc*-antibonding  molecular-orbital  components 
in  Fig.  7  promote  the  covalent  tt  overlap  between  A1  p  orbitals  in  each 
"waveguide  channel,"  thereby  isolating  each  "channel"  from  the  others  much 
like  the  Independent  optical  fibers  in  a  fiber  bundle. 

The  nature  of  the  supercurrent  associated  with  such  a  wavefunction 

8  9 

likewise  follows  from  the  original  work  of  London  and  Slater.  Classically, 
the  electric  current  density  produced  by  n  electrons  per  unit  volume  of 
charge  -e  moving  with  velocity  v  is  given  by  the  formula 

j  *  -nev.  (2) 

In  an  applied  magnetic  field  of  vector  potential  A,  the  electron  velocity 
is  related  to  the  electron  momentum  p  by  the  expression 

m v  =  p  -  eA/c.  (3) 

Combining  equations  (2)  and  (3)  and  reformulating  the  resulting  expression 
in  quantum-mechanical  terms,  the  electric  current  density  associated  with 
a  state  represented  by  a  molecular-orbital  wavefunction  *  occupied  by  n 

g 

electrons  is 

j  =  (-neih/2m)(ij>*V'i'-W*)  -  (ne^mcjA^’V  (4) 

For  example,  one  may  substitute  for  tp  in  the  above  expression  a  simple 
1 inear-combination-of-atomic-orbitals  (LCAO)  representation  of  the  extended 
real-space  superconducting-state  molecular-orbital  wavefunction,  <jj$,  of 
aluminum  shown  in  Fig.  7.  Then  the  terms  with  the  gradients,  which  corres¬ 
pond  to  the  average  value  of  ordinary  electron-momentum  carrying  current, 
cancel  each  other  because  of  symmetry,  leaving  for  the  net  current  density 
the  purely  diamagnetic  term 
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j$  *  -ne2^Us|2/mc.  (5) 

Since  n  is  the  number  of  electrons  in  the  superconducting-state  molecular 

o 

orbitals  of  probability  density  \  ,  the  product 

n$  =  n  1  \ps  | 2  (6) 

represents  the  number  of  electrons  per  unit  volume  which  are  responsible 

for  the  current  density  (5).  Note  that  in  the  cluster  molecular-orbital 

representation  of  a  superconductor,  ip  corresponds  to  only  one  discrete 

energy  level  at  the  Fermi  energy  [the  4t^u(eF)  level  of  the  Al^g  cluster  in 

Fig.  3],  so  that  ng  is  only  a  small  fraction  of  the  total  number  of  electrons 

of  the  system  (the  electrons  in  all  the  occupied  cluster  molecular  orbitals 

8 

of  Fig.  3).  Indeed,  London  has  shown  that,  for  a  simply  connected  isolated 
superconductor  and  for  given  macroscopic  boundary  conditions,  a  single 
electronic  state  of  wide  spatial  extent  (such  as  the  molecular  orbital 
in  Fig.  7)  is  sufficient  to  represent  all  currents  realizable  according  to 
macroscopic  electrodynamics.  Substitution  of  Eq.  (6)  into  Eq.  (5)  yields 

Js  *  -nse2-S/mc,  (7) 

Taking  the  curl  of  both  sides  of  Eq.  (7)  yields  the  magnetic  field  expression 

$  *  v  x  t  =  -(mc/nse2)v  x  js,  (8) 

Q 

which  is  known  as  London's  equation.  London  has  shown  from  Eqs.  (7)  and  (8) 
that  the  current  j$  produced  by  the  magnetic  field  in  a  material  of  macroscopic 
size  is  confined  to  a  thin  surface  layer  and  is  parallel  to  the  surface,  it 
acts  to  shield  the  inner  superconductor  from  the  external  magnetic  field  (the 
Meissner  effect),  and  it  is  not  subject  to  dissipation  or  electrical  resistance 
through  interaction  with  the  lattice.  Thus  the  "supercurrent"  j$  is  of  a 
fundamentally  different  character  from  ordinary  electrical  conduction. 

For  known  superconductors ,  where  Tc  <  25°K,  the  superconducting-state 
molecular  orbitals  <j>s  are  exclusively  occupied  at  the  Fermi  energy  only  for 
T  'v  Q°K.  As  T  approaches  T  ,  thermally  Induced  electronic  excitations 

v 
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gradually  begin  to  depopulate  these  orbitals  and  to  occupy  otherwise  empty 

orbitals  lying  close  in  energy  to  ij»  .  Most  of  the  available  orbitals  above 

the  Fermi  energy  are  antibonding  or  nonbondlnq  in  the  local  chemical  sense, 

as  illustrated  in  Fig.  2  by  the  and  orbitals,  respectively,  of 

the  benzene  molecule.  They  lack  the  coherent  spatially  extended  bonding 

character  of  the  superconducting-state  orbitals  and  therefore  can  be  viewed 

as  constituting  the  nonsuperconducting  phase  observed  for  T  >  Tc.  Since  the 

entropy  of  a  phase  depends  essentially  on  the  number  of  stationary  states 

belonging  to  the  phase,  then  the  entropy  of  the  superconducting  phase,  which 

is  determined  by  a  few  orbitals  i|»  ,  is  negligible  compared  with  the  entropy 

of  the  nonsuperconducting  phase,  which  is  determined  by  the  manifold  of 

excited-state  orbitals.  The  effectively  zero  entropy  of  the  superconducting 

phase  is  consistent  with  the  high  degree  of  order  of  the  superconducting  state, 
g 

Slater  has  argued  that  the  gradual  onset  of  electronic  transitions  between 
the  superconducting  state  and  low-lying  excited  states,  when  translated  into 
the  language  of  thermodynamics,  corresponds  to  the  onset  of  the  second-order 
transition  between  the  superconducting  and  nonsuperconducting  phases. 

Since  the  SCF-Xa-SW  method  satisfies  Fermi  statistics,  the  above- 
described  excitations  can  be  represented  qualitatively  by  the  cluster  molecular 
orbital  model  in  conjunction  with  the  transition-state  concept. 16,17  For 
example,  the  first  allowed  electronic  transition  in  an  Allg  cluster  repre¬ 
senting  superconducting  aluminum  (see  Fig.  3)  is  between  the  pir-po*  molecular 
orbital  4t|u(ep)  shown  in  Figs.  5  and  6  and  the  virtual  dpo  molecular  orbital 
2t2g.  The  2t^g  molecular  orbital  lacks  the  coherent  spatially  extended 
character  of  the  4t^uUp)  orbital  and  is  unoccupied  in  the  ground  state  of 
the  Al-jg  cluster  because  there  are  no  occupied  d  orbitals  in  the  ground  state 
of  an  A1  atom.  However,  this  molecular  orbital  becomes  occupied  even  in  the 
ground  state  when  a  transition-metal  impurity  such  as  Mn  is  substituted  for 
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the  central  atom  in  an  Al^g  cluster,  resulting  in  o-type  chemical  bonding 

between  the  d  orbitals  of  and  p  orbitals  of  nearest-neighbor  A1  atoms,  as 

indicated  in  the  t2g(ep)  wave^unction  contour  map  of  Fig.  9.  It  will  be 

shown  below  that  such  bonding  is  responsible  for  the  degrading  effect  of 

transition-metal  impurities  on  the  superconducting  state  of  aluminum.  For 

pure  Alig,  the  transition-state  energy  of  the  4tluUp)  -*■  2t2g  molecular- 

orbital  excitation  (see  Fig.  3),  which  is  the  discrete  cluster  analogue  of 

the  onset  of  the  phase  transition  between  the  superconducting  and  nonsuperconducting 

states  of  aluminum  particles  and  crystals,  is  ^  0.05  eV,  implying  a  Tc  well 

above  room  temperature  for  this  small  cluster.  Of  course,  it  is  no  more 

practical  to  think  in  terms  of  pir-electron  superconductivity  in  a  19-atom 

aluminum  cluster  than  it  is  in  an  isolated  benzene  molecule,  since  no 

reasonably  sized  coherence  length  is  attained.  Nevertheless,  as  the  cluster 

size  is  systematically  increased  by  including  successive  shells  of  atoms  in 

the  crystal  (e.g.,  in  the  43-atom  results  described  in  Ref.  34),  the  energy 

gap  between  the  tju(£p)  pir-po*  molecular  orbital  at  the  Fermi  energy  and 

first  unoccupied  orbital  decreases  while  the  number  of  unoccupied  orbitals 

immediately  above  the  Fermi  energy  increases,  implying  a  decreasing  Tc  with 

increasing  cluster  size.  Superconductivity  has,  in  fact,  been  observed  for 
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aluminum  particles  down  to  50  A  in  diameter,’  which  is  roughly  an  order 
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of  magnitude  less  than  Slater's  estimate  of  the  coherence  length  (137 

O 

atomic  diameters  =  388  A)  and  many  orders  of  magnitude  less  than  the  coherence 
length  for  crystalline  aluminum  (16,000  A)  derived  from  BCS  theory.  If 
one  extrapolates  the  decreasing  value  of  the  molecular-orbital  energy  gap 

o 

at  the  Fermi  level  with  increasing  cluster  size  to  a  50  A  particle,  using 

a  scaling  procedure  for  the  cluster  energy  levels  analogous  to  that  described 
9  -4 

by  Slater,  a  gap  of  ^  10  eV  is  obtained,  which  implies  a  Tc  of  -v  1°K. 


The  measured  value  of  Tc  for  small  superconducting  aluminum  particles  is 
1.5  to  2.5°K,®  while  Tc  for  bulk  crystalline  aluminum  is  1.14°K.^  As  for 
the  8CS  formula  (1)  for  T  expressed  in  terms  of  the  electronic  density  of 
states  at  the  Fermi  energy,  N(ep),  it  should  be  emphasized  that  only  the 
pir-bonding  molecular-orbital  contribution  to  N(ej-)  (see  Section  II  and  Fig.  4) 
is  pertinent  to  the  superconducting  property  of  aluminum. 

In  the  above-described  chemical  picture,  the  well  known  mutually 
exclusive  incidence  of  superconductivity  and  magnetism  at  a  particular 

3 

temperature  (or  pressure)  in  a  material  can  be  understood  directly  and 
simply  from  the  exactly  opposite  molecular-orbital  criteria  for  the 
occurrence  of  these  two  phenomena.  For  superconductivity,  there  must  be 
a  spin-paired  molecular  orbital  at  the  Fermi  energy  with  coherent  spatially 
extended  bonding  components,  e.g.,  the  pn  components  of  aluminum  described 
above  and  the  ds  components  of  transition  metals  to  be  discussed  below.  In 
contrast,  for  magnetism  there  must  be  a  spin-polarized  spatially  localized 
exclusively  antibondinq  molecular  orbital  at  the  Fermi  energy,  e.g.,  the 
<fa*-dir*  orbitals  of  ferromagnetic  iron  discussed  in  Section  I  (see  Fig.  1) 
and  the  do*-d7r*  orbitals  of  locally  magnetic  Fe  and  Mn  impurities  in  a 
crystalline  copper  host  above  the  Kondo  temperature  discussed  in  Ref.  35. 

The  above  theoretical  model  for  superconducting  aluminum  is  also 
consistent  with  the  observation  that  certain  types  of  magnetic  impurities  in 
aluminum  destroy  its  superconductivity,  or  at  least  lower  Tc.  A  good  example 
is  Mn,  whose  effects  in  aluminum  can  be  modeled  in  a  fashion  analogous  to 
the  treatment  of  Mi  in  copper  (see  Ref.  35)  by  substituting  a  Mn  atom  for 
the  central  A1  atom  in  the  19-atom  cluster  model  of  crystalline  aluminum. 

This  yields  a  MnAl^Alg  cluster  representing  the  local  crystal  environment 
of  the  impurity  up  to  second-nearest  neighbors.  The  Mn  d  orbitals  overlap 
and  hybridize  with  the  surrounding  A1  p  orbitals  of  the  otherwise  unoccupied 
2t2g  orbital  of  pure  aluminum  in  Fig.  3,  resulting  in  the  occupied,  spatially 
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local  1  zed  tggUp)  molecular  orbital  shown  in  Fig.  9.  The  net  overlap  of 
A1  pit  bonding  atomic  orbitals  in  the  t-|u(£p)  orbital,  argued  above  to  be 
essential  for  superconductivity,  is  decreased  around  the  Mn  impurity  due  to 
the  absence  of  localized  p  orbitals  on  Mn,  as  shown  in  the  contour  map  of 
Fig.  10  [cf.  Fig.  6],  and  there  is  a  competition  between  the  t2g(ep)  and 
tlu^eF^  or*)'ita^s  for  electrons.  These  effects  tend  to  degrade  the  super¬ 
conducting  state  by  perturbing  its  long-range  oherercy,  and  the  degrading 
effect  will  increase  with  the  concentration  of  impurities.  Mi  impurities  in 

aluminum,  while  not  producing  permanent  magnetic  moments  below  the  Kondo 
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temperature,  TK  ~  900°K,  give  rise  to  localized  spin  fluctuations.  The 
competition  of  the  degenerate  localized  t2gUp)  dpo  and  extended  t-j u(ep)  P* 
molecular  orbitals  for  electrons  at  the  Fermi  energy,  leading  to  fluctuations 
in  orbital  occupancy,  implies  a  competition  between  the  tendency  to  form 
local  magnetic  moments  and  the  tendency  to  form  the  coherent  spatially 
extended  molecular  state  necessary  for  superconductivity. 

It  is  also  possible  for  molecular  orbitals  of  the  unperturbed  coherent 
bonding  type  essential  for  superconductivity  and  orbitals  of  the  localized 
anti  bonding  type  associated  with  local  magnetic  moments  or  collective 
magnetism  to  be  practically  degenerate  at  the  Fermi  energy.  This  could 
lead  to  a  competition  between  superconductivity  and  magnetism  as  a  sensitive 
function  of  orbital  occupancy,  which  would  depend  on  temperature  according 
to  Fermi  statistics  and/or  on  pressure  according  to  orbital  ordering.  Such 
a  situation  might  explain  the  rare  Incidence  or  "coincidence"  of  super¬ 
conductivity  and  magnetism  in  the  same  material,  e.g.,  chromium  which  is 
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anti  ferromagnetic  at  ordinary  pressure  and  superconducting  at  high  pressure 
(see  discussion  below)  and  superconducting  molybdenum  sulfides  RExMOgSg 
(Chevrel  compounds)  containing  a  high  concentration  (x  ^  1)  of  magnetic 
rare-earth  (RE)  ions.  A  preliminary  molecular-orbital  analysis  of  the 
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MogSg  octahedral  cluster  which  is  the  key  structural  unit  in  the  superconducting 

Chevrel  phases  has  been  carried  out,  using  the  results  of  recent  SCF-Xa-SW 

53 

studies  of  similar  Mo-atom  cluster  compounds  by  Cotton  and  Stanley.  This 

analysis  suggests  that  the  superconducting  state  of  these  Chevrel  phases  and 

54 

the  more  common  ones  not  containing  RE  ions  arises  exclusively  from  coherent 

spatially  delocalized  d-orbital  bonding  between  the  Mo  atoms,  promoted  in  part 

by  the  direct  covalent  overlap  between  Mo  d  orbitals  in  the  MOgSg  clusters 

and  in  part  by  the  antibonding  contribution  of  the  interaction  between  Mo  d 

orbitals  and  S  ligand  p  orbitals.  The  latter  contribution  may  be  compared 

with  the  above-described  effect  of  the  pa *-anti bonding  component  of  aluminum 

in  promoting  coherent  spatially  extended  pit  bonding  (see  Figs.  5  through  7). 

On  the  other  hand,  the  magnetic  properties  of  the  RExMogSg  compounds  are 

associated  exclusively  with  localized  f  orbitals  on  the  RE  ions  which  are 

formally  but  weakly  antibonding  with  respect  to  the  MOgSg  clusters.  Thus  the 

Interaction  between  the  magnetic  RE  ions  and  the  superconducting  electrons 

is  weak,  explaining  why  the  presence  of  such  ions  does  not  generally  destroy 

the  superconductivity  of  the  Chevrel  phases.  More  detailed  SCF-Xa-SW  molecular- 

orbital  studies  of  these  systems  are  planned. 

As  a  somewhat  simpler  illustration  of  the  way  the  d  orbitals  of  transition 

metals  combine  to  form  a  coherent  spatially  extended  superconducting  state, 
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consider  the  A15  compounds,  of  which  NbgSn  is  the  prime  example.  In  this 
compound,  the  Nb  atoms  form  a  network  of  chains  throughout  the  crystal,  with 
the  Sn  atoms  located  at  the  comers  and  center  of  the  unit  cell.  SCF-Xa-SW 
cluster  molecular-orbital  studies  reveal  the  presence  of  ds-bonding  molecular 
orbitals  along  the  Nb  chains  at  the  Fermi  energy.  The  d6-bonding  configura¬ 
tion  is  shown  for  four  atoms  of  a  chain  in  Fig.  11  and  is  consistent  with 
the  type  of  coherent  wavefunction  postulated  above  to  be  essential  for  the 
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superconducting  state.  A  perspective  drawing  of  d<5  bonding  along  a  chain  of 

atoms  is  shown  in  Fig.  12.  Like  the  preceding  example  of  the  Chevrel  phases, 

the  d6  bonding  along  the  Nb  chains  is  promoted  partially  by  the  direct  covalent 

overlap  between  Nb  d  orbitals  (note  that  the  Nb-Nb  interatomic  distance  in 

NbjSn  is  less  than  that  in  elemental  crystalline  nobium)  and  partially  by  the 

interaction  between  Sn  and  Nb  atoms.  The  influence  of  the  Sn  atoms  on  Nb-Nb 

bonding  is  analogous  to  that  of  the  Cl  ligands  on  Mo-Mo  bonding  in  the  molecular 

4+ 

transition-metal  coordination  complex  Mo2Clg  ,  where  the  Mo-Mo  bond  distance 
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is  likewise  shorter  than  that  in  crystalline  metallic  molybdenum.  This  type 
of  ligand-metal  interaction  is  a  significant  contributing  factor  to  the 
variation  of  T  ,  which  is  proportional  to  metal-metal  6-bond  overlap,  among 
the  superconducting  A15  compounds.  It  has  been  customary  to  estimate  T_  values 
for  such  compounds  using  the  BCS  formula  (1)  and  the  total  density  of  states 
at  the  Fermi  energy,  N(ep),  derived  from  band-structure  calculations.57  The 
occupied  d  bands  of  the  A15  compounds  can  generally  be  partitioned  into  do, 
dir ,  and  ds  bonding  components  in  order  of  increasing  energy  toward  the  Fermi 
energy,  with  d6*,  du*,  and  do*  antibondirg  components  lying  above  the  Fermi 
energy.  However,  the  do  and  dn  components  do  not  correspond  to  coherent 
spatially  extended  wavefunctions  of  the  type  argued  above  to  be  necessary  for 
superconductivity  and  therefore  should  be  completely  neglected  when 
calculating  N(ep)  to  estimate  T  . 

The  d  bands  of  elemental  transition  metals  can  also  be  reconstructed 
in  terms  of  overlapping  o-,  n-,  6-bonding  and  6*-,  n*-,  a*-antibonding 
cluster  molecular-orbital  components  and  the  position  of  Fermi  energy  with 
respect  to  these  components  used  as  a  simple  gauge  of  whether  the  metal  is 
likely  to  be  a  superconductor,  magnet,  or  neither  of  these.  For  example, 
the  transition  metals  in  the  vanadium  column  of  the  periodic  table,  which 
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correspond  to  a  valence-electron-to-atom  ratio  e/a  =  5,  are  superconducting 
with  relatively  high  transition  temperatures,  V(Tc  =  5.38°K),  Nb(Tc  =  9.50°K), 
and  Ta(Tc  3  4.48°K),  because  the  Fermi  energy  intersects  the  lower  half  of  the 
d  band  where  there  are  significant  bonding  (e . g. ,  d6)  components  having  the 
coherent  spatial  character  necessary  for  superconductivity.  In  other  words, 
for  these  elements  the  d<$  orbital  contribution  to  the  density  of  states  at 
the  Fermi  energy  is  large.  On  the  other  hand,  although  the  second-  and  third- 
row  transition  metals  in  the  chromium  column  of  the  periodic  table  (corres¬ 
ponding  to  e/a  3  6)  are  also  superconductors,  their  transition  temperatures 
are  relatively  small.  Mo  (Tc  =  0.92°K)  and  W  (Tc  3  0.012°K),  because  the 
Fermi  energy  intersects  the  d  band  near  its  center  where  the  dfi  orbital 
density  of  states  is  small.  This  type  of  analysis  explains  the  variation 
of  Tc  for  transition  metals  (and  their  alloys)  with  electron-to-atom  ratio 

3 

first  pointed  out  by  Matthias.  Chromium  itself  is  normally  anti  ferromagnetic 

50 

and  becomes  superconducting  only  at  higher  than  normal  pressure.  The 
incidence  of  both  anti  ferromagnetism  and  superconductivity  in  Cr  can  be 
explained  by  an  argument  somewhat  analogous  to  that  presented  above  for 
the  coexistence  of  magnetic  rare-earth  ions  and  superconductivity  in  the 
Chevrel  phases.  If  the  Fermi  energy  passes  through  the  center  of  the  Cr 
d  band  where  spatially  extended  bonding  orbitals  of  the  type  necessary  for 
superconductivity  and  spatially  localized  antibonding  orbitals  of  the  kind 
associated  with  magnetism,  i.e.,  the  bonding  and  antibonding  d-band  "edges" 
respectively,  are  close  in  energy,  then  the  orbital  occupancy  and  dominant 
electronic  state  can  be  a  sensitive  function  of  interatomic  distance.  Spin- 
unrestricted  SCF-Xa-SW  cluster  molecular-orbital  studies  of  Cr  for  an 
interatomic  distance  corresponding  to  normal  pressure  do  indeed  indicate 
spin-polarized  anti  bonding  d  orbitals  at  the  Fermi  energy  and  spins  of 
alternating  sign  or  neighboring  atoms,  consistent  with  the  anti  ferromagnetic 
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state.  Similar  studies  of  Cr  for  interatomic  distances  corresponding  to 
high  pressures  are  in  progress. 

Toward  the  right  of  the  periodic  table,  the  ferromagnetism  of  the  first- 

row  transition  metals  (Fe,  Co,  and  Ni)  is  associated  with  the  Fermi  energy 

coinciding  with  spin  polarized  cb*-  and  dir*-antibonding  molecular  orbitals 

33  37 

of  the  type  shown  in  Fig.  1.  *  On  the  other  hand,  preliminary  theoretical 

studies  suggest  that  the  second-  and  third-row  transition  metals,  Ru  and  Os, 
in  the  Fe  column  of  the  periodic  table  are  superconducting,  rather  than 
ferromagnetic,  because  their  hexagonal  close-packed  local  molecular  coordina¬ 
tion  is  associated  with  coherent  spatially  extended  bonding  d  orbitals, 
instead  of  incoherent  spatially  localized  antibonding  d  orbitals,  at  the 
Fermi  energy.  The  second-row  transition  metal,  Pd,  in  the  Ni  column  of  the 
periodic  table  is  neither  superconducting  nor  ferromagnetic  in  its  pure 

crystalline  form  but  is  strongly  exchange  enhanced  by  magnetic  spin  fluctua- 
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tions.  SCF-Xa-SW  molecular-orbital  studies  of  19-atom  palladium  clusters 

representing  the  local  crystalline  environment  up  to  second-nearest  neighbors 

indicate  the  presence  of  degenerate  spin-polarized  ddn*  and  pda*  antibonding 

orbitals  at  the  Fermi  energy.^  Fluctuations  of  electron  occupancy  between 

these  two  types  of  orbitals  are  the  discrete  cluster  analogues  of  the  observed 

magnetic  spin  fluctuations.  SCF-Xa-SW  cluster  molecular-orbital  studies  of 

the  interstitial  alloy,  palladium  hydride  (PdH),61  and  palladium  films 
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containing  radiation-induced  defects,  both  of  which  are  superconducting 
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with  Tc  ^  3°K  or  greater,  are  in  progress.  Preliminary  results  suggest  that 
the  interstitial  hydrogen  and  defects  in  palladium  alter  the  spatial  nature 
of  the  molecular  orbitals  at  the  Fermi  energy,  thereby  quenching  the  spin 
fluctuations  and  inducing  the  superconducting  state. 

The  above  types  of  theoretical  studies  can  also  be  used  to  explain  why 
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some  metals  in  the  periodic  table  are  neither  superconducting  nor  magnetic. 

Among  these,  the  noble  metals  Cu,  Ag,  and  Au  have  been  the  subject  of  much 

discussion  concerning  the  possibility  of  observing  superconducti vity  at 
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extremely  low  temperatures.  SCF-Xa-SW  studies  of  a  19-atom  copper  cluster 

35 

representing  bulk  crystalline  copper  suggest  that  superconductivity  is 
unlikely  to  occur  at  any  temperature  (or  reasonable  pressure).  The  Fermi 
energy  coincides  with  a  molecular  orbital  of  a^g  "spherical"  symmetry  which 
is  primarily  of  the  so*-antibonding  type  and  which  therefore  lacks  the  coherent 
extended  spatial  character  necessary  for  superconductivity  (see  Fig.  4  of 
Ref.  35  for  a  picture  of  this  orbital).  While  there  is  some  antibonding 
d-orbital  contribution  to  this  orbital  due  to  the  hybridization  effect  of 
the  filled  d  band  below  the  Fermi  energy,  it  is  too  small  and  diffuse  to 
give  rise  to  localized  or  collective  magnetic  moments  in  pure  copper.  The 
only  orbital  that  is  capable,  in  principle,  of  forming  a  superconducting 
state  in  copper  is  an  unoccupied  tlu  pir-pcr*  orbital  of  the  type  shown  above 
(Fig.  7)  to  be  occupied  at  the  Fermi  energy  in  aluminum  and  responsible 
for  its  superconducting  state.  Thus  the  question  of  superconductivity  in 
copper  depends  on  whether  this  type  of  orbital,  which  lies  well  above  the 
Fermi  energy,  can  be  brought  into  coincidence  with  the  Fermi  energy  by  any 
practically  attainable  pressure.  Similar  arguments  are  applicable  to  silver 
and  gold. 

While  emohasis  in  this  paper,  thus  far,  has  been  placed  on  materials 

where  the  superconducting  state  is  determined  mainly  by  pir  or  d<$  bonding 

components,  other  types  of  chemical  bonding,  although  less  common,  can  also 

lead  to  coherent  spatially  extended  wavefunctions  of  the  kind  necessary  for 

superconductivity.  For  example,  SCF-Xa-SW  cluster  molecular-orbital  models 
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for  superconducting  "layered-structure"  transition-metal  compounds  such 
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as  NbSg.  TaSg.  NbSe2,  and  TaSe2  indicate  that  the  Fermi  energy  is  coincident 

with  molecular  orbitals  of  dz2(a]'*)  symmetry,  shown  schematically  in  Fig.  13 

for  three  neighboring  metal  atoms  in  a  single  layer.®6  This  orbital  is 

moderately  anti bonding  between  the  transition -metal  d  orbitals  and 

chalcogen  "ligand"  p  orbitals  (not  shown  in  Fig.  13)  but  weakly  bonding 

between  the  "doughnut-like"  annular  regions  of  the  metal  dz2  orbitals  in  each 

layer  (the  dashed  profiles  shown  in  Fig.  13),  in  agreement  with  the  chemical- 
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bonding  model  of  Krebs.  The  bonding  between  annular  regions  forms  a  diffuse 
but  spatially  coherent  wavefunction  throughout  each  layer,  thereby  satisfying 
the  orbital  criterion  for  the  existence  of  superconductivity.  In  this  respect, 
the  superconducting  state  of  the  layered  transition-metal  di chalcogen Ides  is 
"two-dimensional."  However,  the  total  superconducting  state  of  the  crystal 
is  a  composite  result  of  the  three-dimensional  array  of  layered  dz2  orbitals, 
just  as  the  superconducting  states  of  aluminum  and  A15  compounds  discussed 
above  are  the  composite  results  of  three-dimensional  networks  of  "quasi-one- 
dimensional"  p-ir  and  d6  orbital  configurations.  The  relatively  low  transition 
temperatures  of  these  layered  superconductors,  e.g.,  NbS2(Tc  =  6°K)  and 
TaSe2  (Tc  =  0.2°K),  as  compared  with  those  of  the  d6-bonded  A15  compounds, 
e.g.,  Nb3$n  (Tc  =  18°K)  and  VjSn  (Tc  =  3.8°K),  are  consistent  with  the  fact 
that  the  bond  overlap  between  the  annular  parts  o-f  the  dz2  orbitals  in  the 
former  compounds  (see  Fig.  13)  is  significantly  less  than  that  between  the 
d6  orbitals  in  the  latter  compounds  (see  Fig.  12).  The  exact  amounts  of 
bond  overlap  and  concomitant  Tc  values  for  the  layered  superconductors  are 
dependent  on  the  anti  bonding  contributions  of  the  chalcogen  ligand  p  orbitals 
to  the  d^Uj'*)  molecular  orbitals,  in  much  the  same  fashion  as  the  pu 
bonding  orbital  components  in  aluminum  are  promoted  by  the  pa*  antibonding 
components  (see  Figs.  5  through  7),  and  in  much  the  same  way  as  the  net  d6 
bonding  along  the  chains  of  metal  atoms  in  the  A15  compounds,  NbjSn  and  V^Sn, 
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is  dependent  on  the  interaction  of  the  metal  atoms  with  the  neighboring  Sn 

ligands.  The  intercalation  of  organic  molecules,  such  as  pyridine,  between 

the  layers  of  superconducting  transition-ratal  dichalcogenides  does  not 
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appreciably  alter  their  Tc  values.  This  is  explained  by  the  fact  that  the 
bonding  between  the  layers  and  intercalate  molecules  occurs  via  the  dz2  lobes 
perpendicular  to  each  layer  (the  solid  orbital  contours  in  Fig.  13),  which 
leaves  largely  unaffected  the  bonding  between  annular  regions  within  each 
plane  (the  dashed  orbital  contours  in  Fig.  13)  responsible  for  the  super¬ 
conducting  state.  The  weak  annular  type  of  bonding  between  dz2  orbitals  is 
not  strictly  limited  to  layered  transition-metal  compounds  but  can  also 
contribute,  along  with  the  more  dominant  d6  bonding,  to  the  superconducting 
state  in  transition  metals  and  their  alloys. 

In  principle,  spherically  synmetric  s  orbitals  can  also  overlap  to  form 
coherent  spatially  extended  a-type  bonding  molecular  orbitals  like  that  shown 
schematically  in  Fig.  14  and  thus  can  lead  to  the  superconducting  state. 
However,  it  is  unusual  to  find  materials  where  this  kind  of  s-orbital  bonding 
occurs  at  the  Fermi  energy.  It  is  more  common  for  s a  bonding  orbitals  to 
have  energies  well  below  the  Fermi  energy,  e.g.,  the  la-j^  cluster  molecular 
orbitals  and  r-|(^=o)  Bloch  eigenstates  of  aluminum  (Fig.  3),  copper  (Fig.  1 
of  Ref.  35),  and  lithium  (Fig.  3  of  Ref.  67).  As  discussed  above,  the 
molecular  orbitals  at  the  Fermi  energy  in  copper  at  normal  pressure  are  s a* 
antibonding  orbitals  which  cannot  lead  to  the  superconducting  state,  whereas 
pir  bonding  orbital  components  at  the  Fermi  energy  are  responsible  for  the 
superconducting  state  of  aluminum.  The  absence  of  so  or  pn  bonding  orbitals 
at  the  Fermi  energy  of  lithium®^  explains  the  nonsuperconducting  state  of 
this  and  similar  alkali  metals  at  normal  pressures.  Nevertheless,  preliminary 
SCF-Xa-SW  cluster  molecular-orbital  models  for  PdHx  suggest  that  diffuse  so 
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bond  overlap  of  the  kind  shown  in  Fig.  14  occurs  at  the  Fermi  energy 

between  interstitial  hydrogen  atoms,  under  the  chemical  influence  of 

neighboring  palladium  atoms,  for  large  concentrations  (x  ^  1)  of  hydrogen 
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and  may  play  a  key  role  in  determining  the  superconducting  state  of  PdH. 

Finally,  while  the  coherent  spatially  extended  bonding  of  f  orbitals 
has  not  been  addressed  in  this  paper,  their  possible  contribution  to  the 
superconducting  states  of  transuranium  elements  at  normal  and  high  pressure 
are  currently  under  theoretical  investigation. 

IV.  CONCLUSIONS 

Evidence  has  been  presented  in  this  paper  that  a  material  is  likely  to 

be  a  superconductor  if  the  electronic  states  at  the  Fermi  energy,  defined  by 

the  highest  occupied  molecular  orbitals  calculated  by  the  SCF-Xa-SW  method 

for  clusters  representing  the  local  molecular  environment  in  the  material, 

have  coherent  spatially  extended  (e.g. ,  prr,  d6,  dz2,  or  so)  bonding 

components.  It  has  also  been  shown  that,  while  these  individual  orbital 

components  are  usually  "one-dimensional"  or  "two-dimensional"  in  nature 

(see  Figs.  8,  12,  13,  and  14),  the  composite  superconducting  state  of  the 

material  arises  from  the  coherent  three-dimensional  "network"  or  "array"  of 

these  components.  Moreover,  it  has  been  shown  that  for  each  such  orbital 

component,  covalent  bond  overlap  in  a  line  or  plane  of  atoms  is  promoted 

by  the  antibonding  repulsion  of  atoms  not  in  that  line  or  plane,  as 

exemplified  by  the  pir-pa*  interaction  in  aluminum  (see  Figs.  5  through  7) 

and  by  the  ligand-metal  interactions  in  transition-metal  compounds.  This 

results  in  the  type  of  "rigid"  or  "frozen"  wavefunction  of  wide  spatial 

8  9 

extent  originally  postulated  by  London  and  Slater  to  be  necessary  for 


the  production  of  nondissipative  diamagnetic  currents.  These  orbital  criteria 
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can  be  used  to  explain  simply  and  directly  observed  chemical  trends  of 
superconductors  (e.g.,  trends  of  Tc  with  chemical  bonding),  to  rule  out 
certain  types  of  substances  as  possible  superconductors,  and  to  suggest 
ways  of  systematically  improving  existing  classes  or  synthesizing  new  classes 
of  superconducting  materials. 

For  example,  it  has  been  suggested  that  certain  types  of  organic 
molecules  or  polymers  of  more  or  less  one-dimensional  structure  with  conju¬ 
gated  bonds  might  not  only  be  superconducting  but  could  possibly  persist  in 
this  state  at  room  temperature  and  higher.^  Since  three-dimensionality 
seems  to  be  essential  for  the  stability  or  "rigid"  character  of  the  super¬ 
conducting  state,  superconductivity  of  any  kind  and  most  certainly  high  Tc 
superconductivity  in  truly  one-dimensional  systems  can  be  ruled  out.  Further¬ 
more,  while  the  highest  occupied  molecular  orbitals  in  most  conjugated  organic 
molecules  are  *  orbitals  associated  with  the  carbon  p  electrons,  they  are 
typically  of  the  ir^  and  ir^  symmetries  illustrated  for  benzene  in  Fig.  2. 

Such  orbitals  have  nodes  along  the  bond  direction,  e.g.,  the  phase  of  the 
wavefunction  changes,  and  therefore  lack  the  spatial  coherency  necessary  for 
superconductivity.  These  conditions  are  encountered,  for  instance,  in  the 

much  studied  tetrathiaful valene  (TTF)-tetracyanoquinodimethane  (TCNQ) 
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system.  Although  this  substance  in  crystalline  form  exhibits  a  high  degree 
of  ordinary  electrical  conductivity  along  a  certain  crystallographic  direction 
related  to  the  one-dimensional  stacking  of  the  TTF  and  TCNQ  molecules, 
SCF-Xa-SW  studies  indicate  that  the  *  molecular  orbitals  near  the  Fermi 
energy  have  nodal  structures  which  are  the  antithesis  of  the  coherent  spatial 
character  essential  for  the  existence  of  the  superconducting  state. ^  Thus 
the  possibility  of  superconductivity  in  TTF-TCNQ  is  ruled  out  from  first 
principles. 

On  the  other  hand,  crystalline  polymeric  sulfur  nitride  ( SNx )  is  not 
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only  a  good  ordinary  electrical  conductor  along  the  direction  of  "fiber 
bundles"  formed  by  large  aggregates  of  Individual  polymer  molecules,  but  is 
also  an  anisotropic  superconductor,  albeit  a  low  T  0.25°K)  one.70 
SCF-Xa-SW  studies  reveal  that  the  Fermi  energy  is  coincident  with  coherent 
pit  bonding  molecular  orbitals  spatially  delocalized  along  the  (SN)X  polymeric 
chains.71  Despite  the  effectively  one-dimensional  character  of  these  pu 
orbitals,  they  lead  to  a  reasonably  stable  or  "rigid"  superconducting  state 
for  each  fiber  as  a  result  of  pir  interchain  coupling  in  the  (100)  plane  and 
pa  interchain  coupling  in  the  (T02)  plane,71  These  orbitals  are  therefore 
somewhat  analogous  to  the  pir-pa*  molecular  orbitals  shown  earlier  (see 
Figs.  5  through  7)  to  be  responsible  for  the  superconducting  state  of 
aluminum.  However,  in  (SN)X  the  pit  components  lie  exclusively  along  the 
direction  of  the  fibers  and  thus  lead  to  anisotropic  superconductivity, 
whereas  in  aluminum  there  are  three  equivalent  pu  components  along  the  three 
principal  cubic  crystallographic  directions  which  lead  to  isotropic  super¬ 
conductivity.  The  pir  anisotropy  and  relatively  weak  pit  bond  overlap 
associated  with  the  "zig-zag"  structure  of  the  (SN)X  polymers  are  most  likely 
responsible  for  the  relatively  low  transition  temperature  of  this  material. 

The  element  phosphorus,  which  is  closely  related  to  sulfur  and  nitrogen, 

becomes  metallic  and  superconducting  for  moderate  pressure  where  its  crystal 
72 

structure  is  cubic.  SCF-Xa-SW  cluster  molecular-orbital  studies  of  cubic 
phosphorus  suggest  that  for  an  interatomic  distance  corresponding  to  this 
pressure  the  highest  occupied  molecular  orbitals  have  coherent  pit  bonding 
character  along  chains  of  phosphorus  atoms  which  are  coupled  in  a  three- 
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dimensional  network  of  pir  orbital  components  by  the  interchain  interactions. 
This  type  of  bonding  is  responsible  for  the  superconducting  state  of 
phosphorus  (under  pressure),  just  as  it  is  for  (SN)X  and  aluminum.  The 


relatively  high  transition  temperature  (Tc  ^  6°K)  for  phosphorus  is  due  to 
the  fact  that  the  cubic  structure  under  pressure  promotes  coherent  prr  bond 
overlap.  This  suggests  that  the  transition  temperature  of  ( SN)X  would  be 
raised,  in  principle,  if  it  could  be  crystallized  in  a  structure  of  higher 
three-dimensional  symnetry  which  would  also  promote  pu  bond  overlap. 

3 

The  above  observations  provide  strong  theoretical  support  for  Matthias' 

arguments  that  the  highest  transition  temperature  for  a  particular  type  of 

superconducting  material  will  occur  in  the  structural  modification  of  highest 

symnetry.  These  observations  also  suggest  that  in  the  quest  for  higher  Tc 

superconductors,  it  would  be  more  profitable  in  the  long  run  to  attempt  to 

optimize  systematically,  through  structural  or  compositional  modification 

(e.g.,  alloying  or  compound  formation),  the  individual  chemical  bonding 

(e.g. ,  pir,  dfi,  etc.)  components  of  high-symmetry  materials.  For  example, 
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the  existence  of  high  T  superconducting  Chevrel  phases  discussed  above, 

where  the  basic  structural  units  are  symmetrical  clusters  containing  coherent 

metal-metal  bonds,  and  the  existence  of  small-particle  superconductors  with 

5  6 

higher  transition  temperatures  than  the  corresponding  bulk  materials  ’ 
suggest  that  other  compounds  or  composite  materials  based  on  coherently 
bonded  metal  clusters  could  be  potential  superconductors  of  relatively  high 

V 

One  could  go  on  and  give  many  other  examples  where  the  above-described 
molecular-orbital  criteria  for  superconductivity  can  consistently  explain 
the  chemical  trends  of  known  superconductors,  e.g.,  intermetallic  compounds, 
solid-solution  alloys,  and  amorphous  alloys, as  well  as  suggesting  new  types 
of  superconductors.  However,  such  examples  will  be  the  subjects  of  future 
papers.  It  remains  here  to  comment  on  only  one  more  example  where  there 
has  been  much  speculation  on  the  possibility  of  superconductivity,  namely 
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75-77 

metallic  hydrogen  prepared  at  very  high  pressure.  Preliminary  SCF-Xa-SW 

molecular-orbital  models  suggest  that,  while  the  electronic  structure  of  solid 
hydrogen  for  high-pressure  interatomic  distances  and  crystal  structure  is 
consistent  with  metallic  behavior,  this  substance  is  unlikely  to  be  a  super¬ 
conductor.  The  Fermi  energy  coincides  with  an  so*  antibonding  molecular 
orbital  of  the  type  already  shown  above  to  be  responsible  for  the  nonsupercon¬ 
ducting  state  of  copper.  Thus,  as  in  the  example  of  copper,  the  question  of 
superconductivity  in  metallic  hydrogen  depends  on  whether  an  unoccupied 
spatially  coherent  pit-bonding  orbital,  which  lies  well  above  the  Fermi  energy, 
can  be  brought  into  coincidence  with  the  Fermi  energy  by  any  attainable 
pressure. 
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FIGURE  CAPTIONS 


Fig.  1.  Contour  map  of  the  highest  occupied  "spin-up"  molecular-orbital 

wavefunction  of  t2g  symmetry  for  a  15-atom  iron  cluster  representing 
the  local  molecular  environment,  up  to  second-nearest  neighbors, 
in  bcc  crystalline  a-iron.  The  contours  are  plotted  in  the  (110) 
plane  containing  four  of  the  eight  nearest  neighbors  and  two  of  the 
six  second-nearest  neighbors.  The  solid  and  dashed  contours  represent 
positive  and  negative  values,  respectively,  of  the  wavefunction,  with 
the  nodes  between  solid  and  dashed  contours  indicating  dda*  anti- 
bondinq  character  between  nearest  neighbors  and  ddir*  anti  bonding 
character  between  second-nearest  neighbors.  This  spin-polarized 
t2g  wavefunction  is  the  discrete  cluster  molecular-orbital  analogue 
of  states  at  the  top  of  the  majority-spin  d  band  near  the  Fermi 
energy  in  ferromagnetic  crystalline  a-iron.  See  Refs.  26  and  37  for 
further  details. 

Fig.  2.  Comparison  between  the  real -space  molecular-orbital  and  k-space 
band -structure  representations  of  the  pn  electronic  states  of  a 
benzene  molecule,  based  on  the  one-dimensional  periodicity  of  the 
carbon  atoms  (see  Ref.  46).  The  arrows  denote  the  spin  pairing  and 
degeneracy  of  each  molecular  orbital  and  the  "Fermi  energy" 
separates  the  occupied  states  from  the  unoccupied  ones. 

Fig.  3.  SCF-Xa-SW  molecular-orbi tal  energy  levels  of  a  19-atom  aluminum 

cluster  representing  the  local  molecular  environment,  up  to  second- 
nearest  neighbors,  in  fee  crystalline  aluminum.  The  orbitals  are 
labeled  according  to  the  irreducible  representations  of  the  0^ 
cluster  syrmietry  group.  The  4t-|u(ep)  level  corresponds  to  the 
highest  occupied  cluster  molecular  orbital.  Shown  for  comparison 
are  the  energy-band  profiles  for  aluminum,  plotted  along  principal 
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symmetry  directions  of  the  Brillouin  zone. 

Fig.  4.  Electronic  density  of  states  for  a  19-atom  cluster  representing 

the  local  molecular  environment  in  crystalline  aluminum,  determined 
by  a  Gaussian  broadening  of  the  discrete  molecular-orbital  energy 
levels  shown  in  Fig.  3.  The  arrow  denotes  the  position  of  the  Fermi 
energy.  See  Ref.  34  for  details. 

Fig.  5.  Contour  map  of  the  highest  occupied  molecular-orbital  wavefunction 
of  t-ju  symmetry  for  a  19-atom  cluster  representing  the  local 
molecular  environment  in  crystalline  aluminum.  This  orbital  corres¬ 
ponds  to  the  4tlu(ep)  energy  level  shown  in  Fig.  3.  The  contours 
are  plotted  in  the  (200)  plane  containing  four  of  the  twelve 
nearest  neighbors  and  four  of  the  six  second-nearest  neighbors. 

The  solid  and  dashed  profiles  represent  positive  and  negative  values, 
respectively,  of  the  wavefunction.  Regions  of  net  overlap  between 
contours  of  the  same  sign  are  bonding,  whereas  regions  of  positive 
and  negative  contours  separated  by  nodes  are  anti  bonding.  Note  the 
pir  bonding  character  along  the  "horizontal"  direction  and  the 
complementary  po*  anti  bonding  character  along  the  "vertical" 
direction.  Since  this  molecular  orbital  is  triply  degenerate, 
there  are  two  other  equivalent  wavef unctions  with  pir  and  pa*  components 
directed  along  orthogonal  crystallographic  directions,  resulting  in  an 
isotropic  orbital  charge  distribution. 

Fig.  6.  Contour  map  of  the  cluster  molecular  orbital  shown  in  Fig.  5  but 
plotted  in  the  (110)  plane. 
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Fig.  7.  Real-space  molecular-orbital  representation  of  the  electronic  states 
of  crystalline  aluminum  at  the  Fermi  energy,  derived  from  the 
generalization  of  the  t^uUp)  cluster  molecular-orbital  configuration 
shown  in  Fig.  5  to  include  more  shells  of  atoms.  Note  the  coherent 
spatially  extended  pn  bonding  character  along  the  "horizontal" 
direction  and  the  complementary  pa*  anti bonding  character  along  the 
"vertical"  direction.  As  in  Figs.  5  and  6,  the  triple  degeneracy 
of  this  molecular  orbital  implies  two  equivalent  pir-pa*  wavefunctions 
along  orthogonal  crystal  directions,  thereby  resulting  in  an  iso¬ 
tropic  orbital  charge  distribution.  This  three-dimensional  "network" 
of  pir  and  pa*  orbital  components  results  in  the  type  of  "rigid"  or 
"frozen"  wavefunction  of  wide  spatial  extent  originally  postulated 
by  London  (Ref.  8)  and  Slater  (Ref.  9)  to  be  necessary  for  the  super¬ 
conducting  state. 

Fig.  8.  Perspective  drawing  of  spatially  coherent  pw  bonding  molecular- 
orbital  components  along  a  line  of  atoms.  The  solid  and  dashed 
contours  represent  positive  and  negative  values,  respectively,  of 
the  wavefunction. 

Fig.  9.  Contour  map  of  the  highest  occupied  molecular-orbital  wavefunction 
of  *2g  symmetry  for  a  19-atom  cluster  representing  the  local  molec¬ 
ular  environment  of  an  isolated  substitutional  Mn  impurity  in 
crystalline  aluminum  up  to  second-nearest  neighbors.  The  contours 
are  plotted  in  the  (200)  plane  containing  the  Mn  impurity,  four  of 
the  twelve  nearest-neighbor  A1  atoms,  and  four  of  the  six  second- 
nearest  neighbor  A1  atoms.  The  solid  and  dashed  contours  represent 
positive  and  negative  values,  respectively,  of  the  wavefunction. 

The  a-type  chemical  Interaction  between  the  four-lobe  Mn  d  orbital 
and  four  nearest-neighbor  A1  p  orbitals  is  clearly  evident.  Note 


that  this  tgg  molecular  orbital  is  unoccupied  for  pure  aluminum 
in  its  superconducting  ground  electronic  state  (the  2t^g  level 
lying  above  the  Al^g  cluster  Fermi  energy  in  Fig.  3). 

Fig.  10.  Contour  map  of  the  highest  occupied  molecular-orbital  wavefunction 
of  t^u  symmetry  for  a  19-atom  cluster  representing  the  local 
molecular  environment  of  an  isolated  substitutional  Mn  impurity 
in  crystalline  aluminum  up  to  second-nearest  neighbors.  The 
contours  are  plotted  in  the  same  (110)  plane  as  the  corresponding 
orbital  for  pure  aluminum  in  Fig.  6.  Direct  comparison  of  Fig.  10 
and  Fig.  6  reveals  that  the  absence  of  localized  p  orbitals  on  the 
Mn  atom  at  the  cluster  Fermi  energy  reduces  the  net  overlap  of  A1 
Pit  orbitals,  thereby  locally  perturbing  the  coherency  of  the  pit 
bonding  component  and  degrading  the  superconducting  state. 

Fig.  11.  The  spatially  coherent  dfi  bonding  component  of  the  highest  occupied 
molecular  orbital  for  the  A15  compound,  Nb^Sn,  shown  along  four 
atoms  of  a  Nb  chain.  The  three-dimensional  network  of  such  dfi 
components,  coupled  with  the  Sn(ligand)-Nb(metal)  interactions,  is 
responsible  for  the  superconducting  state  of  NbgSn  and  other  A15 
compounds . 

Fig.  12.  Perspective  drawing  of  spatially  coherent  d6  bonding  molecular- 
orbital  components  along  a  line  of  atoms.  The  solid  and  dashed 
contours  represent  positive  and  negative  values,  respectively,  of 
the  wavefunction. 

Fig.  13.  Perspective  drawing  in  a  plane  of  three  atoms  of  spatially  coherent 
dz2  bonding  molecular-orbital  components  of  the  type  responsible 
for  the  highest  occupied  cluster  molecular  orbitals  of  layered 
transition -metal  dichalcogenides  (e.g.,  NbS2).  The  spatially 
extended  coherent  bond  overlap  of  the  annular  regions  (the  dashed 


contours)  of  the  dz2  orbitals  In  each  layer,  coupled  with  the 
effects  of  Nb( metal )-S( ligand)  Interactions,  Is  responsible  for 
the  superconducting  state  of  such  compounds. 

Fig.  14.  Perspective  drawing  of  spatially  coherent  so  bonding  molecular- 
orbital  components  along  a  line  of  atoms. 
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